home *** CD-ROM | disk | FTP | other *** search
/ SGI Developer Toolbox 6.1 / SGI Developer Toolbox 6.1 - Disc 4.iso / src / exampleCode / MP / kernel_switch / bloc_mpy.f next >
Encoding:
Text File  |  1994-08-02  |  726 b   |  36 lines
  1. c 
  2. c Copyright (C) 1993, Silicon Graphics, Inc.
  3. c All Rights Reserved.
  4. c
  5.     subroutine bloc_mpy(c1,c2,m,n,p,alpha,a,lda,b,ldb,beta,c,ldc,
  6.      +   iflag)
  7. c
  8.     implicit double precision(a-h,o-z)
  9.     character*1 c1,c2
  10.     real*8 alpha,beta
  11.     dimension a(lda,p),b(ldb,n),c(ldc,n)
  12.     parameter(iblock=64,kblock=32,jblock=16)
  13.     integer p
  14. c
  15.  
  16.     do j1=1,n,jblock
  17.       j2=j1+jblock-1
  18.       if(j2.gt.n) j2=n
  19.       do k1=1,p,kblock
  20.         k2=k1+kblock-1
  21.         if(k2.gt.p)k2=p
  22.         do i1=1,m,iblock
  23.            i2=i1+iblock-1
  24.             if(i2.gt.m) i2=m
  25. c         call sgi_dgemm_mm( a(i1,k1),lda,b(k1,j1),ldb,c(i1,j1),ldc,
  26.           call mpy( a(i1,k1),lda,b(k1,j1),ldb,c(i1,j1),ldc,
  27.      $                   (i2-i1+1), (j2-j1+1), (k2-k1+1), alpha)
  28.  
  29.     enddo
  30.     enddo
  31.     enddo
  32.  
  33.  
  34.     return
  35.     end
  36.